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During clearonic film 
development, convential film 
is scanned electronically to 
produce a scries of views of 
the dcveling image. An cariy 
scan reveals the fast developing 
highKght detail, white a fate scan 
reveals slow developing shadow 
detail. After development, the 
scries of views is combined 
into a single image in a process 
called stitching. In the present 
I invention, regression data is 
accumulated during development 
to describe a curve of density 
versus time of development for 



I each pixel. After development, 
this regression data is used to 
I recreate a regression curve of 
dentisity versus development 
lime for each pixel. The time 
at which this curve crosses a 
density known to give optimum 
grain characteristics, called the 
optimum density curve, is used 
to create the brighmess for that 
pixel in the finished stitched 
image. The invention further 
teaches weighting regression data 
I as a function of time and density 
; generally following proximity to optiium 
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PARAMETRIC IMAGE STITCHING 

In electronic film development, conventional film is scamied electronically during 
development to produce a series of views of the developing image. An early scan reveals the fast 
developing highlight detail, while a late scan reveals slow developing shadow detail. After 
development, the series of views is combined into a single image in a process called sdtching. In 
the prior art. stitching cut out the best parts of each view and merged them together. In the present 
invention, regression data is accumulated during development to describe a cun^c of density versus 
time of development for each pixel. After development, this regression data is used to recreate a 
regression curve of density versus development time for each pixel. TT^e time at which this curve 
crosses a density known to give optimum grain characterisrics. called the optimum density curve 
is used to cr^te the brightness for that pixel in the finished stitched image. Tbc invention fimher 
teaches weighting regression data as a foncUon of time and density generally following proximity 
to the optimum density curve. 



BACKGROUND AND PRIOR ART 
Recording an image at different exposures aird later merging the images has been practiced 
smce the advent of photography. A technique known to photographers for ove«x,ming the 
dynamic range limit of fihn is to make two exposures, perhaps one for the clouds and one for the 

20 ^hadowedforeground.andmergethetwousinginanualprintingskillinthedarkroom. Asimilar 
technique is known in astrophotography where multiple exposures reveal different feature of a 
star cluster or nebula. In a rather flashy example. Kodak developed a fihn in the 1950's capable 
of recording the million to one brightness range of a nuclear test by making a color negative film 
wherein the three color layers were substituted with three monochrome layers of widely different 

25 sensitivities, each developine in r ol or de vH o p r r w ith n diff ere nt d> o co lor. A gain. U.c uiai . ual 
skill of a daricroom printer was relied on to merge the images into one. A firrther example can be 
found m radiology where images can be made with different x-ray voltages to reveal detail in both 
soft and hard materials, then merging the images together. Modem color film typicaUy uses three 
emulsion coatings for each color, each of a different speed. The three are merged simply by 

30 putting all ti.ee togeti^er in one film, thereby getting some benefit of a layer optimized for a 
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particular exposure, but mixed with the grain of other layers not optimized for that particular 
exposure. 

It was not until the advent of electronic fihn development, as taught in US Patent 
No. 5,519,510 issued to the present inventor, that there was a need to merge multiple exposure 
5 images using producaon-Ievel speed and automation. In electronic film development, the merging 
of images is called stitdiing. The background of electronic film development in general and the 
prior art methods of stitching are now presented as a basis of understanding the background of the 
present invention. 

Turning to Figure 1, a scone 102, portrayed as perceived through the wide dynamic range of 
10 the human eye, has highlights 104, midtones 106, and shadows 108. with details in all areas. A 
camera 1 10 is used to project the scene onto a fihn inside the camera. The scene is perceived by 
the film to consist of points of Ught, each with an exposure value which may be mapped along an 
exposure axis 112. 

The film is removed fix)m the camera after exposure and placed in a developer. In electronic 
15 fihn development, an electronic camera 120 views the film by nonactinic mfiaredJight during 
developmeat. As seen after a short developmoit time of perhaps one minute, the fihn 122 still has 
a low density for shadows 124 and midtones 126. but may optimally reveal hi^li^ts 128. As 
seen by the inftared camera 120, inverting for the negative of conventional fihn, the shadows and 
midtones 130 appear black, while the highlights 132 are seen more clearly than at any later time 
20 in development 

Doubling development time to two minutes, the midtones 140 have progressed to an optimum 
density while the highlights 142 may already be overdeveloped and the shadows 144 may still be 
too low in density to reveal a clear image. The fihn 146 would appear to have good midtone detail 
148. but the highlights ISO are aheady white, while the shadows 152 are still black. 

25 DoubUng development time again to a total of four minutes, the shadows have now reached 

an optimum density, but the other exposures are overdeveloped such that in image 162 they may 
2^]pear white with little detail. 

For each exposure, there is an optimum density of development to reveal the clearest image. 
Clarity may be defined technicaUy as the best signal to noise ratio, where signal is the incremental 

30 change in density with exposure, and noise is ttie RMS deviation in density across a region that 
has received uniform exposure, by convention scanned with a 24 micron aperture. For example 
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at one minute of development time, the midtones 126 typically have too low a density, or are loo 
dark, to have enough of a signal level to reveal detail through the noise of the film and capture 
system. On the other hand, at four minutes the midtones 164 are 'Vashed out", such that not only 
is their contrast, or image signal strength, too low. but the graininess of an oveideveloped silver 
halide emulsion gives a high noise. There exists a development time in between these extremes, 
two minutes in this example, wherein the midtones 140 have developed to an optimum density that 
yields the best signal to noise raUo, or image clarity, for that particular exposure value. In this 
example, the shadows reach optunum clarity at four minutes of development 160, and the 
highlights reach optimum clarity at one minute of development 128. In general, the optimum 
density will be different for different exposures, as in this example wherein the shadows 1 60 reveal 
best clarity at a lower density than the highlights 128. 

After the final caphire of the image on the fibn at four minutes, electronic film development 
has captured optimum images for shadows, midtones. and highlights albeit at different develop- 
ment times. These optimum images must be combined to form a single image with clarity 
throughout approximating the original scene as seen by the wide dynamic range of the human eye. 
TTie process of combining these different parts of the image is called stitching. TTie prior art 
conceived this in the classic sense of merging multiple fihns in a dariooom by cutting out the 
shadows, the midtones. and the highlights, lightening and darkening each so the boundaries 
between regions aligned, then stitching these mulUpIe images together into one. 

The advantage of electronic film development is now more easily understood. In conven- 
tional development the fihn must be stopped and fixed at a selected development time, such as the 
two minute development time of this example. Tlie detail of the highUghts revealed at one minute 
is lost in total daricness as conventional development proceeds. Likewise, the detail that might 
have been revealed at four minutes never had the chance to be bom in conventional development. 
Electronic film development turns conventional fihn into a -'un i vm a l" fi lm tha t can be us e d at 
a wide range of exposure indexes, including veiy high exposure indexes not currently practical. 

InFigure l.thesectionofthedensitycurvearDundtheoptimallydevelopedshadows 160is 
copied as segment 170. Next, the density curve around the optimally developed midtones 140 is 
raised on a base value, or pedestal 172. and copied next to curve 170 as curve 174. Theheightof 
the pedestal 172 is adjusted so the two curves 170 and 174 align. Similarly, the curve around the 
optimally developed highlights 128 is adjusted and raised on pedestal 176 to produce cun^e 178. 
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The process woiks in theory, but in practice, development nonuniformities across the image and 
other spaUally dq,endent nonUnearities made the curves difficult to match across an entire image 
so that the stitched image usually displayed contour, at the edges of stitching regions. Obviously 
an improved method of sUtching was needed to realize the full benefits of electronic film 
development 

Often in electronic film development thfere are more than three exposures made of the film. 
For example, an area array camera may view the fibn continuously, generating hundreds of 
exposures. In the prior art these needed to be combined into a limited number of images to 

conserve memory during the captureprocess. For example, in Figure2the exposures made at one- 
half and one minute, exposures 202 and 2H respectively, could be added with fimction block 206 
to produce a single short development image 208. Similarly, various exposures at other 
development times could be added to yield a middle development image 210 and a late 

development image212.TTiese images would then be aUgned. cut. and pasted together at fim^^^^ 
block 220 to yield the finished image 222. A problem is immediately seen if the times of capture 
vary, making it necessary to adjust the densities in each of the intemrediate images by known time 
deviations. In the past the adjustments were based on estimations of development speed, and were 
not found to be reliable. In addition, there were difficulties if some of the capture times were 
missing entirely because, perhaps, a non real time operating system did not release computer 
resources exactly when needed 

Electronic film development held the promise of higher speed universal film that would work 
in conventional cameras. ITus higher speed and wider range film would enable famiUes to record 
their fives beautifidly in the natural fight of real fife, without typical problems caused by contrast 
bght or dependence on a cold and harsh flash. However, the prior art implem^itations of 
electromc fibn development were plagued with problems in stitching the multiple exposure images 
together. Obviously, an improved stitching meUiod is an important advance to the art 

OBJECTS OF THE IPfVENTION 
The primary object of the invention is to merge images of differing densities into a single 
image which is fi-ee from the artifacts encountered in the prior art. 

A related object is to merge images of differing densities free of edge contouring. 



noise. 
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A further object is to mer^e images of differing densities with reduced effect tor. nonimage 



A further object is to merge images of differing densities while compensating for a shift in a 
density-affecting parameter, such as time. 

to be merged. 

SUMMARY OFTHE INVENTION 
In *= pr=«„, m,««io„. a series of a«cap>™d el«tomcally . d«el.pi„. flta 
usged wi* .he of cap.u:e. p., ^ „f ^^^^^ „ ^ ' 

P=r3„,««s„e«l=„,.,^.s„eh=sde„si.y.tacs toe. or density u„« toe s,„,«>. Ue« 
pamne^. for each faeares^^ed into p.ra:„«„.ecun,u,s,,„g.„,,s. As.refinemem d,. 
pan>m.ttn,canb.»dEhMprior«,s»™u,^hyafl««„rse„s,av.toa,e„Uabai.yof«^ 
FoHomng Sta developmen,. to regression sufedcs are no, neoessari,, vie>«ble i„,ages, ra«,» 
*ey describe ta abstract mathemaUcal ,en„s smooth continuous .i„« for each pixel .ha, pass 

to„gh«,eac««lsan,p.edde„sidesforeachpixelT»ese™ad.«naUcal.ydescribodsmo^^ 
altow d,e dev.,o^„, ^ ^ ^ ^ 

wh,ch U,e densio- of each pixel is p,edic,«. » have attained iB opto™, d««i,y. A gan™, 

oor««ionfa,«io„ofdus„on,„3n««dfaefor.achpixelis.henou,p„,.stobHght„essfor.ha, 
pixel. 

DESCRIPTION OFTHE DRAWINGS 
Figure I depicts the prior art of electronic film development with stitching. 
Figure 2 further portrays prior art sUtching. 
Figure 3 portrays density versus time of a typical development cycle. 
Figure 4 introduces the graphic basis of the invention. 

Figure 5 adds to Figure 4 the effect of measure ment n oise and timebase^ 



Figure 6 adds to Figure 5 the effect of lost data. 
Figure 7 portrays an unweighted embodiment of the invention schematically. 
Figure 8 illustrates a problem with an unweighted embodiment. 
Figure 9 portrays weighting proportional to proximity to an optimum density curve 
Figure 10 presents the preferred embodiment as a series of steps. 
Figure 1 1 presents the preferred embodiment schematically. 
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DETAILED DESCRIPTION OF THE INVENTION 
Sometimes a different way of looking at a problem makes it easier to recognize a new 
solution. The graphs of Figure I described above in the background section followed the prior art 
photographic convention of plotting density as a function of exposure for a series of specific 
5 development times. Figure 3 follows a much less common approach of plotting density versus 
development time for a series of exposures. .Other than that, the points 324, 326, 328, 344, 340, 
342, 360, 364, and 366 in Figure 3 intersect exactly the same triplet of density, development time, 
and exposure as their similarly numbered counteqjarts 124, 126, 128. 144, 140, 142, 160, 164, and 
166 of Figure 1. 

10 The optimum density point for highlights 328, for midtones 340, and for shadows 360 lie on 
a locus of points called the optimum density curve 370. shown by a dotted line in Figure 3. This 
curve is found empirically by measuring signal to noise ratio for vaiying exposures, and finding 
the density at which each reveals detail with the optimum clarity. 

An image can be thought of as consisting of an array of points, or pixels, each of which 
receive a specific exposure. The prior art of Figure 1 thought of this exposure aiH>n>ducing a 
specific density. Further, each development time produced images of differing density that could 
be cut, aligned in density, and merged together. Figure 4 suggests thinking of the image as 
consisting of an array of pixels, each having received a specific exposure resulting in a specific 
development curve, such as the highlight curve 402. Each specific development curve can be 
20 quantified by the time at which the density of the development curve crosses the optimum density 
curve 404. In this case a specific highUght pixel produces curve 402. which crosses curve 404 at 
point 406, which can be tagged as a one minute pixel. 

Continuing with Figure 4, assume a pixel in a middle shadow developed to densities at one. 
two, and four minutes shown by solid dots 410. 412. and 4 14. AlUiough none of these densities 
^ faJIs directly on tiie optimum density curve 404. this way of seeing the problem makes it clear tiiat 
a best fit line 418 could be drawn through the known points 410 to 414 to predict that Uie pixel 
crossed tiie optimum density curve 404 at point 416. so that pixel can be tagged as a 2.8 minute 
pixel. This method requires no cutting and aligning of multiple images; rather the placement of 
the best fit line provides a continuum between areas of differing exposure, completely eliminating 
30 edge contouring and thereby solving a significant problem in the prior art. 



15 
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Further, assume that there was a time base error in sampling the data. In the prior art of 
Figure 2. an attempt to coirect the error would have involved estimating what effect that error had 
on the early, middle, and late images. In the case of Figure 4. however, the error would appear as 
a shift in the sample times, for example 1.5. 2.5. and 4.5 minutes shown by the x's 420, 422. and 
424. It may be noted that the shifted sample times would still lie along line 418. and therefore the 
best fit line to points 420. 422 and 424 and the best fit line to 410. 412 and 414 is the same line 
418 which still crosses the optimum density ctirve at the same point 416 at 2.8 minutes. In fact, 
individual sample points could be added, deleted, or moved in time with minimal effect on the best 
fit curve or the estimate of the crossover time, thereby solving another significant problem in the 
prior art. 

Figure 5 illustrates a more typical case wherein the electronic camera adds noise to the 
captured images. In this case, each capture will not only lie along a Une that is a fimction of film 
exposure and grain, but in addition each capture adds a random noise deviation fiom the line. TTie 
individual density samples for a particular pixel are illustrated in Figure 5 as solid dots, such as 
dot 502 measured at a time of 2.3 minutes. The density of sample 502 differs fiom a theorized tme 
curve 504 because of noise in the electronic camera arising from, typically, statistical errors in 
counting photons, called shot noise. Combined with several other samples, however, illustrated 
in Figure 5 as the multiple solid dots, a best fit curve 504 can be estimated, and the time this best 
fit curve 504 crosses the optimum density curve 506 calculated, as before. 

A particular question arises fiom Figure 5 as to how to specify a best fit curve, arid the answer 
providesanotherdistinctioninthepracticeoftheinventionoverthepriorart. Well known statistic 
practices give an array of choices. One option is to gather regression parameters on the points, and 
then calculate a linear best fit fiom Uie sum of densities and ttie sum of densities times time: By 
including the sum of densities times time squared as one of the parameters, a quadratic regression 
can be used to yield a quadratic fimcUon for t h r, best fit curve Re^t fit regre s sion analysis of such 
parameters is well known in the art. For example, a linear best fit is a curve of the forai density= 
A+Bt, where A and B are found such that the sum of the square of Uie distances of each sampled 
point Smm the line is minimized. 

Now one of Uie distinctions in the practice of the invention over the prior art can be stated. 
Unlike the prior art which sought to gather actual images from the film at specific zones of 
development time and later merge those real images, the present invention seeks to gather more 
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abstract regression parameters about the developing image and later use those parameters to 
recreate a real image for any conceivable development time. The mathematical expression of the 
image as a function of a time continuum then allows a seamless stitching. 

A particularly interesUng best fit function is the locus of lines representing the developed 
density versus time as a function of exposure level for similar film, such as curves 510. 512. 514 
and 516. The best fit curve could be described mathematically as density=F(t.tc). where t is any 
development time at which density is to be solved, and tc specifies the shape of the curve to cross 
the optimum density curve 506 when t=tc. The curve with the best fit to the sampled data, curve 
504 in this example, is selected, and the crossover time read directly as the variable. 

The best fit curve could also be described mathematicaUy as F(t,e). where e is the exposure 
level yielding a particular curve. This locus of lines can be derived by actually developing a series 
of test fihns given known exposures, and storing the actual measured densities of each as a 
fimction of development time t in a lookup table where one axis of the lookup table is the known 
exposure, a second axis is the time since developer induction that a specific measureraait is made, 
and the value stored in the lookup table is the empirically measured density. To distinguish such 
a curve repertoire fiom a repertoire of mathematically simple curves, such as a series of lines in 
linear regression or cun.es described by a quadratic fomiula in a quadratic regression, such a 
regression will be called an empirical curve regression after the locus of curves derived from 
empirically measuring film during development Any specific curve can be read fiom the lookup 
table by selecting a particular exposure value as one axis of the lookup table, and then varying 
time, the other axis, while reading out density. The value of exposure yielding the cunre with the 
best fit is the best fit curve. 

The parameters gathered to specify one of the repertoire of available curves in the lookup 
table can be derived by summing the density of each sample poin t times fiinctions of time A fi,.. 
parameter is the product of a first fimction of time, a second parameter is a second fimction of 
time, and so forth. For example, in a quadratic regression, one of the fimctions of time is time 
squared. These fimctions of time are derivable fit>m the empirical curv«.. The first fimction can 
be related to an average of the curves. The second fimction can be the primary mode by which the 
curves differ fiom this average, called the residue after the first fimction has been subtracted from 
each curve. The third fimction can be the main remaining residue after removing the components 
of the first and second fimctions. and so forth to whatever order is desired. 
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Anolher interesUng parameter set is based on a gaussian function of time, where the gaussian 
function is taken to be a function of time that rises and falls smoothly in a bell shape to select a 
particular period of time. A series of parameters based on such overlapping gaussian functions 
would specify the stnooth shape of a cuive. The use of such curves may be found in the art. 
especially in spatial transforms where the human retinal neural system is found to respond to 
gaussian and difference of gaussian (DOG) functions. 

Turning now to the next illustration. Figure 6 illustrates the case of a real time operating 
system that terminated capture prior to reading all the data. Even though data, as represented by 
solid dots, was not received up to the optimum density time, nevertheless the best fit method 
allows a curve to be found and projected through the optimum density curve. Although this would 
not produce the most grain fiee image, it would produce an acceptable image under conditions in 
which the prior art would have straggled because there were no middle and late exposures to stitch. 
This ability to recover fix>m a failure is again a significant advance over the known art. 

Now that the basis of the invention is understood, a specific embodiment using quadratic 
regression is presented schematically in Figure 7. Figure 7 inputs the same sequence of seven 
scamied images 702 and requires the intermediate storage of three accumulating arrays 704. 706. 
and 708. as shown in the prior art example of Figure 2. The distinction over the prior art is that 
in Figure 7. tiie accmnulating arrays sum regression statistics rather than images. This is 
emphasized by portraying tire three accumulating arrays 704. 706 and 708 with crosshatching to 
indicate they are not meant necessarily as viewable images, but rather as statistical.data. 

In Figure 7 a series of images 702 is received sequentially fiom an electronic camera viewing 
the developing film. For example, image 720 is received at two minutes of development The 
density of each pixel of image 720 is summed with the density of corresponding pixels of images 
taken at other development times, such as image 722 at five minutes. This summation occurs in 
fimction hlnr .k 724 w hi ch can ei t h er o perat e on a ll ima ge s toge l l.er if they have all been 
accumulated and stored in memory during development, or one by one as they are captured from 
the electronic camera at the corresponding development times. The advantage of summing them 
as they are captured is that less memory is required. The resuhing sum from fimction block 724 
is stored in the accumulating array 704. 

Continuing the process of the present invention, the density of each pixel of image 720 is 
multiplied by time, two minutes in this example, and the product summed with the density times 
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time of corresponding pixels of images representing other development times, such as image 722, 
for which t=5 minutes. This summation occurs in function block 730. and the summaUon from 
fimcUon block 730 is stored in the accumulating 706. Finally, the products of density times 
time are multiplied by time again to yield density times time squared, and the summation derived 
at function block 732 is stored in die accumulating array 708. In tiie alternative, tiie process could 
continue to a cube or higher orders of time to support a cubic or higher order regression, or could 
end before Uie square term to support a linear regression. 

After the last image is captared and its statistics summed to tire accumulating arrays, then a 
best fit curve is derived for each pixel by retrieving the corresponding statistical data for thlt pixel 
from tiie accumulating arrays 704. 706 and 708. The time of intersection of U,e best fit curve for 
each pixel witii the optimum density curve is calculated in function block 740, and a function of 
this time stpred for tiie corresponding pixel in tiie fmished image anay 742. TTie fimction stored 
in tiie final image array 742 can be tiie time directiy; or it can be tiie exposure known empirically 
to develop to tiie particular time, found empirically and expressed in tiie computer as a function, 
such as tim,ugh a lookup table, of tiie time; or it can be any otiier function found te-have utility, 
such as tiie square root of tfie linear exposure value nomialized and stretched to fit white level and 
black level, yielding a conventional 8 bit computer image representation. 

The direct parametric embodiment just presented does have limitations fliat are now 
highlighted witii reference to Figure 8. In this figure, a specific highlight exposure curve 802 is 
plotted as density versus development time. A series of noisy samples are represented by solid 
dots such as dot 804 representing a noisy measurement of density at one minute, and dot 806. 
representing a noisy measurement of density at five minutes. The goal is to find flie crossover time 
808. one minute in tins example, where curve 802 crosses Uie optimum density curve 810, witii 
th e added constraint tiiat this point 808 be found witi^out exacflv knowing curve 807. nnW fU. 
noisy measurements such as 804 and 806. 

An immediate problem which is evident is Uiat a linear regression curve, such as straight line 
812, will be significantly affected by samples distant fiom tiie optimum density time, such as 
sample 806. and as a result the best fit linear curve 812 will intersect tfie optimum density curve 
810 at an incorrect time 814. If tiiis merely added a bias enx>r. it would not be a significant 
problem; however, changes in tfie outlying samples, such as 806, can be seen to change tiie cross- 
over time 8 14. and because outlying samples such as 806 bear almost no real information, this 
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change is an unwanted noise. It should be noted that higher o«ler regressions, such as the 
quadratic of Figure 7. would put less emphasis on outlying points; however, higher order regres- 
sions require more accumulating arrays than a linear regression and correspondingly mor. memory 
and computing power, and the reduction in emphasis of outlying points is not precise or complete. 

A solution to this problem is presented in Figure 9. Basically, the statistics for each sample 
are throttled by a weighting factor. In this example, sample 904 (equivalent to sample 804 of 
Figure 8) is given a high weight, while sample 906 (equivalent to sample 806 of Figure 8) is given 
a very low weight. In this specific example, weighting factors are assigned such that a linear 
regression yields straight line 912. which is much closer to the tme curve 902 in the region around 
the crossover time 922. at the expense of deviating at outlying times, such as six minutes for 
which accuracy is irrelevant n,us. one object of this invention to reduce noise is accomplished 
via application of a weighting factor. 

The weighting function should generally foUow the reliability, or signal to noise ratio, of each 
sample. The crest of the weighting fimction should therefore follow the optimum density curve 
shown in Figure 9 as the 100% weighting curve 930. and fall offin proportion to distance fiom 
this curve. To illustrate the application of this weighting process, sample points such as 904 that 
are close to the optimum density curve are shown large, while sample points such as 906 that : 
distant are drawn faintly. In this specific example, only the 100% and 50% weighting curves : 
shown, although nonnally, the feUoff would be a continuous function of distance. 

The rate of falloff with distance follows the rate of falloff of overall signal to noise ratio 
mcluding both film grain and noise in the electronic capture system, and is found empirically. If 
the capture system is very noisy, the feUoff is slower so more points are averaged. If the capture 
system is noise fiee or the film has a narrow range of optimum densities, then the weighting falloff 
IS rapid. In practice there is a wide tolerance, and a curve such as shown in Figure 9 will be very 
close to optimum for a rang ^ nf films. 

Figure 9 fiulher portrays Uie true curve 940 for a shadow exposure with the actual sample 
pomts from which U,e true curve is to be matirematically estimated. Again, sample points such 
as 942 that are proximate to the optimum density curve 930 are portrayed large to indicate a high 
weighting, while outlying sample points, such as 944. are drawn faintly. Again it may be 
visualized that a best fit linear line through the large points around cur^e 940 will closely 
approxunate curve 940 around five minutes at which time the curve crosses the optimum density 
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poinL Curve 940 is included to point out that a weighting factor proportional to distance fiom the 
optimum density curve will tend to emphasize samples at short development Umes with highlight 
exposures, and samples with long development times for shadow exposures. 

Finally, it should be noted that the weighting factor should not drop to zero at very low 
densities or short times, but should maintain some fmite weight. This insures that even if the 
system fails to capture samples at later -firaes, there will be some regression data in the 
accumulating arrays from which a reasonable guess at the best fit curve can be extracted. 
Therefore, sample point 944 would be given a small but non-zero weight. The weight can. 
however, go to zero for high densities and long times. 

Figure 10 portrays Uie embodiment just described as a series of steps such as would be 
implemented in a computer program. The specific example given is one of linear regression. 

Initially, ttiere are available to tiie process tiuee arrays to hold the regression statistics. All 
the elements of tiiese regression arrays arc initially set to zero. Three functions are also available: 
a first function receives a density and time as parameters and returns a weight factor, a second 
fimction receives time as a variable and returns an optimum density, and a third function receives 

acrossover time asaparameter and retuinsabrightness value artistically representativeoftheU^^ 
the crossover time represents, which is typically Oie square root of linear brightness, called a 
gamma correction. 

During development, an electronic camera views tiie film. The image fiom the camera is 
di^tized into pixels. As Uie process continues, each point fiom the film is sensed at several times. 
The process receives these pixels one at a time, receiving for each pixel four numbeis: a density. 
X and y coordinates of Uie pixel, and a time at which tiie specific density was measured. 

Based on these four numbers for each pixel, regression parameters are calculated and 
co llected. Specifically, regression parameters are calculated based on density and time, and array 
elements in the regression summation arrays pointed to by x and y are incremented by tiie 
regression parameters. 

Afier tiie last pixel is received for tfie last time, the summation arrays are complete, and the 
data tiiey hold can be used. For each x, y location, regression data is read fiom ttie arrays and used 
to calculate a linear best fit curve. The time at which tius best fit curve crosses tiie optimum 
density curve is calculated, or read fiom a lookup table, or solved by iteration. Finally. Uie equiva- 
lent x,y element of a final image is set according to a gamma fimction of this crossover time. 
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Figure 1 1 covers the procedure of Figure 10 schematically. A series of images are received 
from an electronic camera viewing a film at specific times since a development has been induced 
For example, image 1 102 is received for two minutes of development. A fiinction such as 11 04 
labeled 'Veighf receives the time and also receives the numeric density of each pixel in the 
image, returning for each pixel a first value that sums to the corresponding pixel in a "summation 
weight" array 1106 via conduit 1108 terminating, in 1110. This first value is also sent to 
multiplication block 1112 along with the numeric density of each pixel to produce an output 
product "summaUon density" for each pixel. "Summation density" sums to the corresponding 

pixel in "summaUondensit/'arraylll4viaconduitni6andni8.AnothermultiplicationbI^^ 
1 120 also receives this "summation density" for each pixel, further receives the corr^ponding 
time, and outputs a product called "summation time" for each pixel that is the product of weight, 
density, and time. "Summation time" sums to the corresponding pixel in the "summation time- 
array 11 22 via conduit 1 1 24 and 1 126. lire process is repeated for each new image, such as image 
1 128. that is received by the electronic camera captured at a different development time. 

Following development, the three parametric arrays 1106. II14. and 1122 are used to 
calculate a brightness value for each pixel. The three values received for each pixel ftom the three 
arrays can be used to regress to a linear equation of the fomi density=A+Bt. and t allowed to vary 
iteratively to solve for the time at which a fimction optimum_density(t)=A-fBt. A faster method 
divides the two values from the summation weight anay and summation time array by tiie 
corresponding value fiom the summation density array, and the resulting two numerical values are 
used as a pointer into a two dimensional lookup table that holds Uie precalculated time, or a 
gamma function of the time. Finally, the value so found is placed into U.e corresponding pixel of 
the finished image 1 128. which now. with tire fiiU brightness perception range of Ute human eye. 
is seen to be Shiriey and little Albert frolicking in tiie bluebonnets. 

^ ' '^ > he prese nt invention ha s b tc u dc^ci ibul liereia with reference to particular 

embodiments thereof, a latitude of modification, various changes and substitutions are intended 
in tire foregoing disclosure, and it will be appreciated tirat in some instances some features of the 
invention will be employed wiUiout a corresponding use of oti,er features without departing from 
the spirit and scope of the invention. 
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CLAIMS: 

1 A method of monitoring change with time of an image comprising a plurality of pixels 
possessing densities that are functions of time, comprising the steps of: 
receiving a first density of a pixel at a first time, 

accumulaUng to a fust parameter sum a first fiincUon of the first density and the fust time, 
accumulating to a second parameter sum a second function of the fust density and the first 
time, 

receiving a second density of the pixel at a second time, 

accumulaUng to the first parameter sum the first fimction of the second density and the second 
time; and 

accumulating to the second parameter sum the second fimction of the second density and the 
second time. 

2. The method of claim I, wherein the first fimction has a nonzero value at the fiisttime and the 
second time. 

3. The method of claim 1, wherein the first fimction is a product of a density and a function of 
time. 

4. The method of claim 3. wherein the fimction of time is time squared. 

5. The method of claim 3, wherein the fimction of time is time. 



6. The method of claim 3, wherein the fimction of time is derived from differences in empirical 
curves of density versus development 

7. The method of claim 3, wherein the fimction of time is a gaussian fimction of time. 

8. The method of claim 3. wherein the first fimction further includes weighting by a fimction 

of the density. 
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9. The method of claim 8. wherein the weighting is a function of density and time. 

10. The method of claim 9. wherein for a selection of values the function has greater weightings 
for values closer to an optimum density curve. 

11. The method of claim 8. further comprising the step of accumulating the weighting to a third 

parameter sum. 

12. The method of claim 1. wherein the plurality of pixels is associated with an associated 
plurality of first parameter sums and an associated plurality of second parameter sums. 

13. The method of claim 12. further comprising the step ofderiving a best fit line to a continuum 
of densities of a pixel as a function of time, wherein the deriving step uses the firet and the second 
parameter sums associated with the pixel. 

14. The method of claim 13. wherein the deriving step includes a polynomial regression. 

15. The method of claim 14. wherein the deriving step includes a linear regression. 

16. The method of claim 14. wherein the deriving step includes a quadratic regression. 

17. The method of claim 13. wherein the deriving step includes an empirical regression. 

18. The method of claim H , w heceinJhe d e riving ste p in cludes a gaussian regr c ^iou. 



19. The method of claim 13. further comprising the steps of: 
selecting an optimum density curve, and 

solving for a crossover time at which the best fit line crosses the optimum density curve. 
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20. The method of claim 19, further comprising the step of creating an image in which a value 
of the pixel is a function of the crossover time. 

2 1 . A method of monitoring change with time of an image comprising a plurality of pixels, 
each pixel having a density that is a function of time, comprising the steps of: 

receiving densities of the pixels at a plurality of times; 

accumulating to a Orst parameter sum for each pixel a first function of the density and the 
time; and 

accumulating to a second parameter sum for each pixel a second function of the density 
and the time. 

22. The method of claim 21, v^^erein the plurality of times includes three distinct tunes. 

23. The method of claim 21, wherein the first function is a product of a density and a function of 
time. 

24. The method of claim 23, whereui the first function further includes weighting by a fimction 
of the density. 

25. The method of claim 24. wherein the weighting is a fimction of density and time. 

26. The method of claim 25. wherein for a selection of values the fimcUon has greater weightings 
for values closer to an optimum density curve. 



27, The method of claim 24, fiirther comprising the step of accumulating the weighting to a third 
parameter sum for each pixel. 



28. The method of claim 21, further comprising the step of deriving a best fit line to a continuum 
of densities for each pixel as a fimction of time, wherein the deriving step uses the first and the 
second parameter sums associated with each pixel. 
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29. The method of claim 28. wherein the deriving step includes a polynomial regression. 

30. The method of claim 28. further comprising the steps of: 
selecting aq optimum density curve, and 

solving for a crossover time at which the best fit line for each pixel crosses the optimum 
density curve. 

3 1. -me method of claim 30. further comprising the step of creaUng an image in which a value of 
each pixel is a function of the crossover time. 

32. The method of claim 21 . forther comprising the step of accumulating to a third parameter sum 
for each pixel a third function of the density and the time. 

33. The method of claim 32. further comprising the step of deriving a best fit line to a continuum 
of densiti^ for each pixel as a fimction of time, wherein the deriving step uses the first, the second 
and the third parameter sums associated with each pixel. 



34. The method of claim 33. wherein the deriving step includes a polynomial 



regression. 



35. The method of claim 33, further comprising the steps of: 
selecting an optimum density curve, and 

solving for a crossover time at which the best fit line for each pixel crosses the optimum 
density curve. 

36. The method of claim 35. fiirther comprising the step of creadng an image in which a value of 
-each^ixel4s- a function of th e cto ssoverttme: 
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